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Abstract 

Under the direction of NASA's Office of 
Aeronautics and Space Technology (OAST), the NASA 
Lewis Research Center has initiated an in-house 
thermal energy storage program to identify com- 
binations of phase change thermal energy storage 
media for use with a Brayton and Stirling advanced 
solar dynamic (ASD) space power systems operating 
between 1070 and 1400 K. A study has been initi- 
ated to determine suitable combinations of thermal 
energy storage (TES) phase change materials (PCM) 
that result in the smallest and lightest weight 
ASD power system possible. To date the heats of 
fusion of several fluoride salt mixtures with 
melting points greater than 1025 K have been veri- 
fied experimentally. The study has indicated that 
these salt systems produce large ASD systems 
because of their inherent low thermal conductivity 
and low density. It is desirable to have PCMs with 
high densities and high thermal conductivities. 
Therefore, alternate phase change materials based 
on metallic alloy systems are also being consid- 
ered as possible (TES) candidates for future ASD 
space power systems. 

Introduction 

Solar dynamic systems are actively being 
studied as a means of providing heat and electri- 
cal power for space applications. Advanced Solar 
Dynamic (ASD) power systems offer the potential 
for efficient, lightweight, relatively compact 
and long-lived space power systems applicable to 
a wide range of power levels (3 to 300 kWe), and 
a wide variety of orbits. Thermal energy storage 
(TES) will be required by these ASD power systems 
to provide power to the heat engine during solar 
eclipse. The latent heat TES concept that will 
be discussed in this paper results from a solid 
to liquid phase change which offers a lightweight 
and direct method of providing the heat engine 
with a uniform heat source. While many critical 
technology issues exist for ASD power systems, 
finding an efficient means of energy storage is 
presently the least well defined and understood. 
Under the direction of NASA's Office of Aeronau- 
tics and Space Technology (OAST), NASA Lewis has 
initiated an advanced solar dynamic testing and 
development program. The goals of this in-house 
program are to identify and select high temper- 
ature (1025 to 1400 K) phase change thermal energy 
storage materials which will yield light weight, 
compact and efficient solar dynamic Brayton and 
Stirling space power systems with an operational 
1 ife of 7 to 10 yr . 

Literature surveys and analytical screening 
has been performed to identify PCM systems with 
heats of fusion greater than 400 kJ/kg. Compar- 
isons have been made between salt and metallic 
PCM systems and it was found that the metallic 
PCMs produced the smallest and lightest weight 


ASD systems because of their higher thermal con- 
ductivities and densities. This paper will detail 
the methods used to identify and select candidate 
PCMs for advanced solar dynamic space power sys- 
tems. Also, the effects of PCM selection on a 
35 kWe Closed Brayton Cycle (CBC) and a Free 
Piston Stirling Engine (FPSE) power system will 
be discussed. 

Analytical Screening 

The primary focus of the in house "Advanced 
Solar Dynamic Systems Code," is to determine the 
impact of the thermophysical properties of the 
selected PCMs on the power conversion systems 
(PCS) mass and performance. Although, selection 
of PCMs with high melting points will increase the 
system efficiency of a Stirling or Brayton power 
system, the system mass was found to also be sen- 
sitive to other thermophysical properties of the 
PCM. The computer systems code used in this study 
will give the minimum system mass for a Brayton 
and Stirling power system operating at the melting 
point of various PCMs.* The goal of the systems 
effort is to identify systems with specific 
weights less than 60 kg/kWe for LEO based service. 

Receiver Considerations 

As shown in Fig. 1, the heat receiver contri- 
butes greatly to the mass of a solar dynamic power 
system. A typical solar dynamic heat receiver is 
shown in Fig. 2(a) and a detail of one PCM/working 
fluid tube is shown in Fig. 2(b), The study was 
done using three salt based PCMs (LiF, LiF-CaF 2 , 
NaF) and two metallic PCMs (Mg-Si and Ge) selected 
from the initial literature survey. Figs. 3 and 4 
show the system masses of a 35 kWe Brayton and a 
35 kWe Stirling PCS using both salt and metal PCMs. 
One critical difference in the masses of the salt 
and metal systems is the sizing of the heat 
receiver cavity. It is desirable to have uniform 
flux distribution within the cavity. Therefore, 
the receiver was sized by the diameter and spacing 
of the storage tubes based on a reasonable maximum 
allowable flux. The larger receiver of both cav- 
ities was used in the analysis. 

One characteristic of salt PCMs is low ther- 
mal conductivity while the metal lies are very good 
conductors. Because of the high thermal conduct- 
ivities of the metal PCMs, higher solar fluxes 
can be tolerated within the cavity resulting in a 
smaller heat receiver, see Fig. 5(a). To improve 
the thermal conductivity of the salt system fins 
and reflectors are required to evenly distribute 
the solar flux within the cavity, this produces a 
large and massive heat receiver, see Fig. 5(b). 

The mass of several of the heat receiver com- 
ponents are directly affected by an increase and 
or decrease in the density and heat of fusion of 
the PCM. For example, if the density were 
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decreased more containment material would be 
required and the weight of the support structure 
and the outlet manifold would increase. Moreover, 
if the heat of fusion is increased the amount of 
storage material required will decrease thus 
decreasing the receiver mass. Therefore it may 
be logical to assume that in order to minimize 
the receiver weight, a PCM with a high density 
and high heat of fusion or a combination of PCMs 
with heat of fusion values high enough to offset 
the mass penalty imposed by low density values 
must be selected. 

Systems Computer Code 

The computer code used in this study deter- 
mines for a given peak cycle temperature the opti- 
mum temperature ratio that will give the minimum 
system mass for the Stirling and Brayton power 
systems. Moreover, the analysis for the Brayton 
power system further optimizes to give the pres- 
sure ratio at the minimum mass conditions. 

Several assumptions were made in generating the 
mass models for the Brayton and Stirling (PCS) 
and the PCMs. Namely, steady state performance 
of both of these systems was assumed. It was also 
assumed all of the heat which is absorbed by the 
PCM is transferred to the PCS working fluid. The 
specific mass of the concentrator and heat pipe 
radiator were assumed to be 1.22 and 5 kg/m 2 
respecti vely. The emissivity of the radiator was 
assumed to be 0.88. 

Analytical Results 

To determine the effects of a PCM's heat of 
fusion and density on the mass of a Stirling PCS, 
a parametric study was done. In this study, the 
density and heat of fusion of Li F was used to 
assess the mass penalties associated with varying 
these thermophysical properties over a specified 
temperature range. First, the heat of fusion was 
varied over a range of 400 to 1200 kJ/kg with the 
density being held constant at 1874 kg/m 3 . As 
shown in Fig. 6, the system mass decreases with 
increasing heats of fusion evaluated at a constant 
density. Similarly, the system mass decreases 
with increasing density when evaluated using a 
specified heat of fusion (Fig. 7). In both 
figures, it should be noted that the rate of 
change in system mass is less at the higher heats 
of fusion and densities. 

Figure 8 shows different scenarios of the 
same system evaluated using combinations of PCMs 
with low, medium, and high densities with high and 
low heats of fusion. Using mass as the selection 
criterion, it is evident that the system with the 
low density and low heat of fusion has a high 
mass. Moreover, the system with high heat of 
fusion and density is the least massive. The 
figure also indicates that there are combinations 
of the heats of fusion and the density where the 
system mass is practically the same. Therefore, 
one of the key factors that will limit the ability 
to identify systems with low masses is the selec- 
tion of PCMs with high heats of fusion and densities 
or with the right combinations of thermophysical 
properties such as thermal conductivity. 

Mass savings of a power system can be sub- 
stantial when selecting PCM systems with both high 
density and high heats of fusion. Figure 9 shows 
that significant weight reduction can be achieved 
by using the ASD systems over the state of the art 


(SOA) systems and Space Station Technology. 2 To 
date LiF has been identified as the best salt 
system for achieving a lightweight system because 
of its high heat of fusion. Yet, this study 
indicates that metal systems produce lightweight 
ASD systems because of the higher thermal conduct- 
ivities and high densities. 

Although the preliminary analysis of these 
systems appear promising, the issue of transferr- 
ing heat into and out of the storage material and 
containment of these PCM systems have yet to be 
addressed. 

Thermophysical Property Determination 

As indicated in the previous section, the 
thermophysical properties of interest to this 
program include heat of fusion, thermal conduct- 
ivity and density. To date, literature searches 
and other sources such as, universities and gov- 
ernment labs have indicated that the thermophysical 
properties are well defined for pure salts but 
not for eutectic salts or metallic alloys. 

NASA Lewis has the capability to determine 
the heat of fusion. Table 1 shows the estimated 
heat of fusion for several fluoride salt mixtures 
in the 1025 to 1400 K temperature range. The 
heat of fusion for many of the fluoride salts has 
now been experimentally verified by differential 
scanning calorimetry (DSC). Table 2 shows the 
heat of fusion of several metallic alloys also 
with melting points between 1025 and 1400 K. 

Many heat of fusion values for these metallic 
alloys have been taken from Barin and Knack. 3 

Other thermophysical properties which must 
be considered to determine the PCMs effect on the 
ASD system are thermal diffusivity and specific 
heat. Currently an in-house capability is being 
developed to determine density and thermal con- 
ductivity, and the capabilities of several labora- 
tories throughout the United States are being 
investigated as possible sources for determining 
the specific heat and or the thermal diffusivity 
of the selected PCMs. 

Technical Issues 

The systems code has identified two types of 
phase change thermal energy storage systems with 
the potential of meeting the lightweight and high 
performance requirements imposed by the ASD power 
system. However, the low thermal conductivity of 
the salt systems and the containment and compat- 
ibility issues associated with both must be 
addressed. 

Thermal Conductivity Enhancement 

As previously stated one inherent character- 
istic of the salt PCMs is low thermal conductivity 
which is currently being enhanced by fins and 
reflectors and results in a large and massive ASD 
power system. One approach that can be used to 
enhance the thermal conductivity is presently 
being investigated at NASA Lewis. In this method, 
if the thermal energy storage canisters of the 
heat receiver shown in Fig. 10 are filled with a 
matrix of high thermal conductivity fibers, such 
as graphite, then significant increases in thermal 
conductivity are possible. The location of the 
fibers with respect to the thermal energy storage 
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salt is critical to both the thermal conductivity 
enhancement as well as the location of the voids 
that form on freezing. If the fibers are wetted 
by the phase change liquid in specific locations 
then enhanced thermal conductivity can be main- 
tained. For example, with the addition of as 
little as 4 % (by volume)of chopped P-100 graphite 
fibers or vapor grown fibers, the thermal conduc- 
tivity of LiF should increase by a factor of 2 to 
6 . 

Containment and Compatibility 

Initial compatibility screening has concluded 
that mild steel, pure nickel and Nb-lZr have 
relatively good corrosion resistance in molten 
fluoride salts . 4 Unfortunately, only Nb-lZr 
possesses some measure of creep strength at ele- 
vated temperatures. However, it is a refractory 
metal alloy which requires special handling and 
fabrication techniques for successful employment. 
There is therefore much reluctance to use this 
type of material. Ideally the strength of the 
containment vessel should be sufficient to with- 
stand the stresses developed during melting of 
the solid PCM. Many of the fluoride salts have 
very large expansions and volume changes upon 
melting for example, LiF increases in volume about 
30 percent when melted. Such behavior could prove 
disastrous if trapped solid is melted, because the 
mechanical forces resulting from the volume change 
could breech the containment vessel wall. Through- 
out its life the heat receiver's TES unit will be 
subject to tens of thousands of freeze-thaw cycles 
under microgravity conditions which pose ill- 
defined solidification and remelting conditions. 

It is desirable for the containment vessel to 
possess adequate strength to cause the solid salt 
to fracture as opposed to itself. For this to be 
the case designers must have knowledge of the 
strength of the PCMs more specifically the salts 
because of the large volume changes they experi- 
ence upon melting. 

Based on the apparent ability of essentially 
pure Fe, Nb and Ni to resist attack, it is pos- 
sible that an alloy containing only these three 
elements would also be resistant to the molten 
salts. This offers the potential for developing 
alloys with inherent high temperature strength, 
for they can be built on the Ni— N i 3 Nb system 
(with various amounts of Fe) that exhibits 
behavior similar to Ni — N i 3 AI which forms the 
basis for nickel superalloys. 

Although, metallic alloys based on Si and Ge 
possess excellent thermophysical properties, they 
react with all metals at high temperature and rule 
out metallic alloys as suitable container materials. 
Ceramic materials such as SiC and Si 3 N 4 are pos- 
sible candidates for containment materials, and as 
indicated by Ref. 5 research is currently being 
conducted to produce a Si containment shell on a 
Si-base eutectic alloy sphere. It is evident to 
all working in the area of containment for Si and 
Ge-based alloys that much developmental work needs 
to be done to find suitable container materials. 

Prior to the final selection of a PCM for the 
ASD power system, long term (5000 to 10 000 hr) 
compatibility testing with potential containment 
materials must be conducted. Thus, the selection 
of a phase change thermal energy storage material 
for ASD power systems is based on analytical and 
metallurgical screening. 


Concluding Remarks 

A study has been initiated to determine how 
PCM thermophysical properties affect the mass and 
size of an ASD power system. This study has shown 
that high temperature, high heat of fusion and 
high density phase change thermal energy storage 
materials must be utilized in order produce small 
and lightweight ASD systems. To date salt and 
metal phase change systems have been investigated. 

As a result of the higher thermal conductivity and 
density of the metallic phase change systems, as 
compared with the salt systems, smaller ASD 
systems were produced. 

A literature survey to determine the thermo- 
physical properties of the desired PCMs has 
indicated that such information is not readily 
available, or that it does not exist. Therefore, 
capabilities are being developed in-house to 
determine the density and thermal conductivity of 
the PCMs and other sources are currently being 
investigated as possible sources for determining 
the specific heat and or thermal diffusivity. 

There are several technical issues associated 
with the PCMs identified as possibly being able to 
meet the lightweight and minimum mass requirements 
of NASA Lewis' ASD program. For the salt PCM the 
major issue is low thermal conductivity and to this 
end an in-house program is under way to investi- 
gate the possibilities of enhancing the thermal 
conductivity with graphite fibers. Another area 
of concern for both the salt and the metal systems 
is containment and in the near future a contain- 
ment development program will be undertaken. For 
the salt systems a series of Ni— Ni 3 Nb alloys con- 
taining varying amounts of iron will be investigated 
and for the Si and Ge-based metallic systems ceramic 
and other new technologies will be explored. 
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TABLE 1. - PHASE CHANGE TES CANDIDATES 
(FLUORIDE SALTS) 1025 TO 1400 K 


Salt 

Melting 

Heat of 

molejt 

temperature. 

fusion. 

K 

kJ/kg 

LiF-22CaF 2 

KF-15MgF 2 

1039 

1043 

a 790 

b 500 

Li 3 A1F 6 

1059 

540 

NaF-32CaF 2 

1083 

590 

NaF-23MgF 2 

1103 

620 

Li F 

1121 

°1087 

NaF-14AlF 3 

1161 

539 

NaF-27CaF 2 -36MgF 2 

1178 

520 

CaF-50MgF 2 

1253 

610 

NaMgF 3 -30KMgF 3 

1259 

641 

NaF 

1268 

b 790 

NaF-60MgF 2 

1273 

670 

KF-68.7MgF 2 

Ha 3 AlF 6 

1281 

1283 

a 690 

b 530 

NaMgF^ 

1303 

c 7 1 0 

MgF 3 

1343 

d 710 

KCaF. 

1343 

c 500 

FeF2^ 

1373 

b 550 

MgF 2 

1536 

b 930 


a Heat of fusion data for MgF£ obtained from 
in-house experimental measurements. 

Stained from Barin and Knacke's 
text .3 

c Heat of fusion data obtained by 
Bukhalova et al J 

^Obtained from in-house experimental 
measurements.® 
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TABLE 2. - PHASE CHANGE TES 
CANDIDATES (METALLIC ALLOYS) 
1025 TO 1400 K 


Composition, 
mole % 

Melting 

temperature, 

K 

Heat of 
fusion, 
kJ/kg 

Ge 

1213 

510 

Mg-54.5Si 

1223 

a 750 

GaAs 

1238 

600 

ZnP 2 

1253 

730 

ZnAs 2 

1288 

450 

Ca-69Si 

1296 

b 1110 

InP 

1328 

430 

Mg 2 Si 

1373 

1118 

Mg 2 Ge 

1388 

705 


a Obtained from Ref. 3. 
b Oak Ridge National Lab data. 



FIGURE 1. - COMPONENT BREAKDOWN OF AN ADVANCED SOLAR 
DYNAMIC POWER SYSTEM. 
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1000 1100 1200 1300 


INLET TEMPERATURE/ K 

FIGURE 3. - SPECIFIC MASS AS A FUNCTION OF TEMPERATURE 
FOR A 35-kWe BRAYTON SOLAR DYNAMIC POWER SYSTEM USING 
VARIOUS PCMs. 
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(A) RECEIVER SIZING BASED ON LiF (SALT) AND Mg-Si (METAL). 


REFLECTORS • 



*- SALT PCM 



-FINS 


GAS FLOW TUBES - 1 L METAL 

PCM 

METAL 

(HIGH CONDUCTIVITY) 
(B) EFFECTS OF LOW AND HIGH THERMAL CONDUCTIVITY. 

FIGURE 5. - COMPARISON OF SALT AND METAL PHASE CHANGE THERMAL 
ENERGY STORAGE SYSTEMS. 


SALT 

(LOW CONDUCTIVITY) 
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HEAT OF 
FUSION, 
KJ/KG 



INLET TEMPERATURE/ K 

FIGURE 6. - SYSTEM MASS AS A FUNCTION OF INLET TEMPERATURE 
FOR A 35-kWe STIRLING SOLAR DYNAMIC POWER SYSTEM USING A 
SALT PCM AT VARIOUS HEATS OF FUSION. NUMBER OF TUBES/ 
60; NUMBER OF FINS/ 6; DENSITY/ 1874 kg/m 3 . 
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INLET TEMPERATURE, K 

FIGURE 7. - SYSTEM MASS AS A FUNCTION OF INLET TEMPERATURE 
FOR A 35-kWe STIRLING SOLAR DYNAMIC POWER SYSTEM USING A 
SALT PCM AT VARIOUS DENSITIES. NUMBER OF TUBES, 60; NUM- 
BER OF FINS, 6; HEAT OF FUSION, 1044 kJ/kg. 
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HEAT OF DENSITY. 
FUSION. KG/M 3 
kJ/kg 



TEMPERATURE. K 


FIGURE 8. - SYSTEM MASS AS A FUNCTION OF INLET TEMPERATURE 
FOR A 35-kWe STIRLING SOLAR DYNAMIC POWER SYSTEM USING 
PCMs WITH VARIOUS HEATS OF FUSION AND DENSITIES. 
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CONCENTRATOR SLOPE ERROR 



INLET TEMPERATURE, K 

FIGURE 9. - SPECIFIC SYSTEM POWER AS A FUNCTION OF INLET 
TEMPERATURE FOR 35-kWe ASD SYSTEM. 




FIGURE 10. - BASELINE SOLAR DYNAMIC HEAT RECEIVER SHOWING THERMAL 
ENERGY STORAGE CONTAINMENT CANISTER FILLED WITH THERMAL- 
CONDUCTIVITY-ENHANCING GRAPHITE FIBERS. 
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